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Abstract: The precise mechanism that causes HIV infection to progress to AIDS is still 
unknown. This paper presents a mathematical model which is able to predict the entire 
trajectory of the HIV/AIDS dynamics, then a possible explanation for this progression is 
examined. A dynamical analysis of this model reveals a set of parameters which may produce 
two real equilibria in the model. One equilibrium is stable and represents those individuals who 
have been living with HIV for at least 7 to 9 years, and do not develop AIDS. The other one is 
unstable and represents those patients who developed AIDS in an average period of 10 years. 
However, further work is needed since the proposed model is sensitive to parameter variations. 
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1. INTRODUCTION 

Several mathematical models have been proposed to de- 
scribe HIV dynamics since 1990, these present a basic 
relation between CD4-I-T cells, infected CD4+T cells and 
virus Nowak [2000], Kirschner [1996], Perelson [1999], Xia 
[2007]. These models give a good presentation of the initial 
peak infection and the asymptomatic stage. However, they 
are not able to describe the transition to the last stage of 
the disease AIDS (acquired immunodeficiency syndrome). 

To obtain a more widely applicable model, some authors 
have tried to introduce other variables, taking into consid- 
eration other mechanisms by which HIV causes depletion 
of CD4-I-T cells. Numerous theories have been proposed, 
but none can fully explain all events observed to occur 
in practice. Recent studies Wang et al. [2000] have shown 
that HIV infection promotes apoptosis in resting CD4-I-T 
cells by the homing process. This mechanism was mod- 
eled in two compartments by Kirschner et al. [2000], in 
this study authors showed that therapeutic approaches 
involving inhibition of viral-induced homing and homing- 
induced apoptosis may prove beneficial for HIV patients. 
The role of the thymus in HIV-1 infection was considered 
by Kirschner et al. [1998]. The authors found that infection 
of the thymus can act as a source of both infectious 
virus and infected CD4-I-T cells. Significant effort has been 
developed in understanding the interaction of the immune 
system and HIV Campello [1999], Adams et al. [2004]. 
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One limitation of these mathematical models is that they 
do not reproduce the entire trajectory of HIV/AIDS dy- 
namics. This trajectory consists of the early peak in the vi- 
ral load; a long asymptomatic period and a final increase in 
viral load with a simultaneous collapse in healthy CD4-I-T 
cell count during which AIDS appears. 

A number of studies have been conducted to explore 
the role of macrophages in HIV infection as long-term 
reservoir Oreinstein [2001]. A reservoir is a long-lived cell, 
which can have viral replication even after many years 
of drug treatment. Using this theory, Conejeros et al. 
[2007] proposed a deterministic model which describes 
the complete HIV/AIDS trajectory. Simulations results for 
that model emphasize the importance of macrophages in 
HIV infection and progression to AIDS, but no dynamical 
analysis is proposed. 

In this paper, we present a simplification of Conejeros et al. 
[2007], in order to have the same behavior of HIV/ AIDS, 
which permits us to understand the transition to AIDS. 
The model is discussed and compared with clinical data. 



2. MODEL DESCRIPTION 

The model proposed in this section is a simplification of 
Conejeros et al. [2007], and considers the following popula- 
tions; T represents the uninfected CD4-I-T cells, Tj repre- 
sents the infected CD4-I- T cells, M represents uninfected 
macrophages. Mi represents the infected macrophages, 
and V represents the HIV population. 



The mechanisms consider for this model are described by 
the next reactions; 

A. Cell proliferation 

The source of new CD4+T ccUs and macrophages from 
thymus, bone marrow, and other cell sources is considered 
constant. 

0^T (1) 
^ M (2) 

si and S2 are the source terms and represent the generation 
rate of new CD4+T cells and macrophages. However, when 
pathogen is detected by the immune system, a signal is sent 
in order to become more aggressive, and then CD4+T cells 
and macrophages proliferate; 

T + V ^ {T + V) + T (3) 

M + V {M + V) + M (4) 

B. Infection cell 

HIV can infect a number of different cells; activated 
CD4+T cells, resting CD4+T ccUs, quiescent CD4+T 
cells, macrophages and dentritic cells. For simplicity, just 
activated CD4+T cells and macrophages are considered 
viral hosts: 

T + V^Ti (5) 

M + V^Mi (6) 

C. Virus proliferation 

The viral proliferation is modeled as occurring in activated 
CD4+T cells and macrophages. 

Ti^V + Ti (7) 
Mi-!^V + Mi (8) 

D. Natural death 

Cells and virons have a finite lifespan. These losses are 
represented by the following reactions; 



T 


-^0 


(9) 


T- 


-^0 


(10) 


M 


A0 


(11) 


Mi 


-^0 


(12) 


V 


A0 


(13) 



Using reactions (1)-(13), we obtain the following model; 

T = si + kiTV - k2TV - 6iT 

ti = hTV - S2T, 

M = S2 + kzMV - kiMV - S3M (14) 

M, = kiMV - SiMi 

V = k^Ti + keMi - 55V 

The model implementation outlined in (14) will be con- 
ducted using MATLAB. Parameters and initial conditions 



wore obtained from previous works in the area Perelson 
[1999], Xia [2007], Conejeros et al. [2007]. Using clinical 
data for the CD4+T cell counts Greenough [2000], Fauci 

ct al. [1996], some parameters were adjusted, see Table 1, 
in order to obtain the best match with clinical data. 

Table 1. Parameters Values 



Parameter 


Value 


Value taken from: 


Sl 


10 


Perelson [1999] 


S2 


0.15 


Perelson [1999] 


ki 


2 X IQ--^ 


Fitted 


k2 


3 X IQ-^ 


Fitted 


ki 


7.45 X IQ-* 


Conejeros et al. [2007] 


fc4 


5.22 X 10"* 


Conejeros et al. [2007] 


fes 


5.37 X 10"^ 


Conejeros et al. [2007] 


fee 


2.85 X 10"^ 


Conejeros et al. [2007] 


Sl 


0.01 


Conejeros et al. [2007] 


S2 


0.44 


Fitted 


S3 


0.0066 


Fitted 


Si 


0.0066 


Fitted 


55 


2.4 


Xia [2007] 



3. IVIODEL RESULTS 

An infected HIV patient, in general, suffers a fast drop 
in healthy CD4+T cell count and at the same time a 
rapid increase in virus population. Then, the immune 
system responds to the virus by proliferating CD4+T 
cells and macrophages, this can be seen just after the dip 
during primary infection in Fig.l. During the next 8 to 
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Fig. 1. CD4+ T cells dynamics. Comparison with clinical 
data taken from Greenough [2000] and Fauci et al. 
[1996] 

10 years the patient experience an asymptomatic phase. 
On one hand, CD4+T cells experience a slow depletion 
but are with sufficient level to maintain most immune 
system functions. At the same time the virus population 
continues infecting healthy cells, therefore slowly advances 
in numbers, see Fig. 2. At the end of asymptomatic period, 
constitutional symptoms appears when CD4-I-T cells are 



below about 300 mm~^. The last stage and the most 
dangerous for the patient is when CD4+T cells drop below 
200 mm~^ and the viral explosion takes place, which is 
considered as AIDS. We can see in Fig.l how the proposed 
model is able to represent all the stages of the infection and 
matches well clinical data. 



Fig. 4 shows how macrophages climb over the years slowly 
trying to suppress the virus. This fact is consistent with 
previous works Conejeros et al. [2007], who suggest that 
macrophages may divide and become more aggressive. 
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Fig. 2. Viral dynamics 

Infected CD4+T cell behavior is structurally similar to 

the viral dynamics in the first years, sec Fig. 3. There is an 
initial peak of infected CD4+T cells, followed by a small 
increment but almost constant during the asymptomatic 
stage. 
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Fig. 3. Infected CD4+ T cell dynamics 

Macrophages play a central role in this model of HIV 
infection. They are considered one of the first points of 
infection and then infected macrophages are long-lived 
virus reservoirs as is noted in Oreinstein [2001]. We may 
notice how the model is able to represent these facts. 
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Fig. 4. Macrophages and infected macrophages dynamics 

Infected macrophages increase slowly in number during the 
asymptomatic period, but when constitutional symptoms 
appear infected macrophages increase very rapidly as can 
be seen in Fig. 4. This is consistent with the work of 
Igarashi et al. [2001] , who argued that in the early infection 
the virus replication rate in macrophages is slower than 
replication rate in CD4+T cells, but over the period of 
years, the viral replication rate in macrophages is faster 
than early stages of infection. 

4. AIDS TRANSITION 

The model proposed by Conejeros et al. [2007] shows the 
complete HIV trajectory. However, they do not give a 

detailed explanation of the HIV/AIDS transition, since the 
model is difficult to analyze. For this reason, we propose 
some simplifying assumptions to allow analysis. 
Assumption.! Fast Viral Dynamics 

Looking at the differential equations (14) and parameter 
values in Table 1, we notice that 5^ >> 1. In this case, the 
differential equation for the virus can be approximated by 
the next algebraic equation, as noted in Barao et al. [2007] . 

y=^T, + ^M, (15) 
Assumption. 2 Assume Ti is bounded 

We note that in the asymptomatic period of infection 
(that is, after the initial transient, and before the final 
divergence associated with development of AIDS) , the con- 
centration of infected CD4+T cells is relatively constant. 
This assumption is also proposed in Astolfi et al. [2008]. 
Therefore the following assumption for infected CD4+T? 
cells can be considered: 



Ti(i)«Ti,Vi>io 
Then under (15), (16) can be reduced to 



(16) 



V{t) := CiMi + Vt, (17) 
where Vt^ = IfTi and ci = f^. Note that if fi is 
selected as an upper bound on Tj, then (17) represents 
an upper bound on v{t). Therefore (17) describes the long 
asymptomatic period in the viral load dynamic. Using last 
assumptions in macrophages and infected macrophages 
equations, we have the following system 

M « S2 - C2M + csMMi (18) 
Mi w C4M + csMMi - S^Mi (19) 
where C2 = ^3 — (fcs — A;4)Vtj, C3 = (fcs — ki)ci, C4 = kiVjr 
and C5 = k4,c\. 

Assumption. 3 M and Mi have an affine relation 

Note that from (18) and (19), we expect that the bilinear 

terms are predominant for large M and Mj, then we may 

assume M w ^Mi, which are rearranged in the next linear 

form; 

Mi W CqM - C7 

where ce = and 07 = ^. 

C3O4 ' ni^A 



(20) 



C3S4 



Proposition 1. Under Assumptions 1-3, the macrophage 
dynamics in an infected HIV patient are unstable with a 
finite time escape. 

Proof: Substituting (20) in equation (18), we have a good 
approximation for the macrophage equation; 



M = S2 + aM^ + i3M 



(21) 



where 

05 



P ~ & ^3 

The solution of the differential equation (21) for Aas2 > 
is given by; 



^ = 2^ + 2^ \ 2 * + ^ J ^''^ 

where 77 is a constant related to initial condition of 
macrophages given by; 



T] = tan' 



(23) 



In (22) there is a tangent function, which tends to 00 when 
the argument tends to w/2, that is when; 



t Trjc .- 



IT — 2r] 



-\/4as2 - P"^ 
implies that there is a finite escape time. 



(24) 



5. DYNAMICAL ANALYSIS 

The structure of the system (14) can be decomposed in two 
feedbacks, see Fig. 5; one is a fast negative feedback and the 
other is a slow positive feedback. The biological meaning of 
this is that feedback 1 is due to a fast CD4+T cell infection 
which is greater than the CD4+T cell proliferation. The 
feedback 2 is for a slow macrophages infection rate which 
is less than macrophage proliferation. 



Feedback 1 




Fig. 5. HIV scheme 

Using the system (14), we are able to get the equilibrium 
points analytically in the form; 



T 



M: 



Si _ fesi V 

kdV + 5i' '~ 62 kdV + Si 

S2 k,S2 V 



knV + ^3 



h k„V + 6s 



where kd = k2 — ki, kn = k^ — k^ and the value of F is a 
solution of the polynomial. 

aV^ + bV'^ + cV = (25) 
The equation (25) has three solutions, which are; 

ViA) ^ y(B) ^ :± + V^Z^^v^c) ^ -b-Vb''-4ac 
2o ' 2a 

where; 

a = 526ii5zknkd 

b = 52^3^4<55fcd + SiS2S465kn — 54^k2knk^Si — 52k4,kdkQS2 

c = 5i5253,5a55 — (53(54^2^531 — kike6i52S2 
Equilibrium A 

T(A) = £1^ t(^) = 0, M(^) = ^, ) = 0, F(^) = 
01 03 



Equilibrium B,C 

j,{B,C) 



Sl 



T- 



(B,C) _ 



S2 



82 fciF(^'t^) + (5i 

(B,C) ^ ^232 V^'''^^ 

5i k2ViB,c)^Ss 



Equilibrium A represents an uninfected status. Using 
numerical values, the uninfected equilibrium is unstable, 
which is consistent with previous works Astolfi et al. 
[2008]. This could explain why it is difficult to revert 
a patient once infected, back to the HIV-free status. 
Using parameter set in Table 1, equilibria B and C take 
imaginary values since 6^ < 4ac, which are not important 



for the biological case. This means that there is no stable 
point for the patient with this set of parameters, then the 
model will progress to AIDS status. 
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Fig. 6. Parameter space in terms of the number of equilib- 
rium points: shaded region shows two equilibrium 

5. 1 Number of equilibrium points and parameter spaces 

Using the recently developed symbolic real algebraic ge- 
ometry methods we keep some parameters unfixed and 
get the number of equilibrium points for the unfixed pa- 
rameter space. The method is called the discriminant va- 
riety method which was developed in Lazard et al. [2007] , 
this is implemented in Maple as in-built packages called 
"Parametric" and "DV" Liang et al. [2009]. This method 
decomposes the parameter space into different cells such 
that the number of real roots is the same for any point in 
a given cell. We compute 2D parameter spaces below, for 
/c4 and k^. In Fig. 6 the colored region shows the part of 
the parameter space in the — k^ plane for which there 
are two equilibria in the system. All other parameters are 
set as in Table 1. 

Remark 1. Using a built-in symbolic algebra routine 
we confirm that the number of equilibria with biological 
meaning and dynamic properties for the proposed model 
are equal to those in Conejeros et al. [2007] . 



5.2 Bifurcation Analysis 

Using the proposed model, there is a set of parameters 
which may produce real equilibria in B and C. If we solve 
for parameter Js, the change of stability is given by; 

^ ^ S2k5kdkQS2 - 5^k2knk^si ^ 
^ {-Sikn + kd5z)525i 

2{-S4k2knk5SiS2k4kdkeS2)^ ^^g-j 

{Sikn + kdS3)S2S4 

Increasing the value of as in (26), the escape time may 
be delayed. This means that infected patients that are 
able to rapidly clear the virus could postpone for many 



years the transition to AIDS. Furthermore, a change in 
the stability from unstable to stable is shown in Fig. 7, 
that is the immune system would adapt to maintain a 
stable status in the patient. 




Fig. 7. Bifurcation using (the black region is the 
unstable behavior and the blue region is stable) 



Table 2. Bifurcation 



Parameter 


Critical Point 


Sl 


50.679 


S2 


2.463 


ki 


2.95x10"^ 


k2 


3.33x10-^ 




5.337x10"^ 


ki 


5.648x10-"* 


k5 


2.721 


k(i 


4.679 


Sl 


1.178x10-^ 


& 


4.732x10^1 


53 


6.99x10"^ 




7.4x10-^ 


S5 


2.497 



Remark 2. Whilst the model reproduces known long 
term behavior, bifurcation analysis in Table 2 evidences 
an unusually high sensitivity. In particular, small relative 
changes in ^2, fcs, k^, Si, 62, S3, or S5 give bifurcation to a 
qualitatively different behavior. 



5.3 Numerical Results 

Choosing ^5 as stated in condition (26), there are two real 
positive equilibria which are shown in Table 3. One point 
is called LTNP (Long Term Non Progressor), because 
it shows the characteristic of individuals who have been 
living with HIV around 7 to 12 years and have stable 



CD4-I-T counts around 600 cells /mnT' or more. The other 
point is called progressor; when CD4-I-T cell counts fall 
below 200 cells /mrn?, the patient is said to have AIDS. 

We observe in Table 4 that LTNP point is stable and the 
progressor point is unstable. The model can represent ade- 
quately both phenomena as is presented in clinical studies. 
Prom the bifurcation analysis, it would appear that there 
are certain progressor patients who could become LNTP. 
However, it is unclear if this property would be preserved 
in more complete models, and clinical evidence for this 
seems weak. 

Table 3. Equilibria 



Points 


V 


T 


T 


M 


M, 


LTNP 


5.43 


647.99 


24.0 


27.83 


11.95 


Progressor 


21.31 


319.3 


46.40 


81.22 


136.92 



Table 4. Eigenvalues for equilibria 



LTNP 


Progressor 


-0.66 


-0.66 X 10-^ 


-0.49 X 10-^ + 0.14 X lO^^i 


0.286 X 10"^ 


-0.49 X 10"^ - 0.14 X 10-^i 


-0.457 X 10"^ 


-0.286 


-0.238 


-3.35 


-3.19 



6. CONCLUSIONS 

The proposed model predicts the entire trajectory of the 
disease: initial viremia, latency, and the rapid increase 
of virus. Using this simplified model can be understood 
how HIV works in an infected patient; basically, the virus 
inhibits the CD4-I-T cell population while promotes the 
macrophages proliferation which are reservoir for virus 
replication. The long reservoir behavior of macrophages in 
an infected HIV patient is a possible explanation of why 
a patient progress to AIDS, and serve as recommendation 
for clinical study. However, we are concerned about the 
sensitivity of the model for some parameters, alternative 
mechanisms will be considered in future work. 
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